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The Angular Impact Distibution of 
Charged Particles Attracted to a 
Charged Cylindrical Spacecraft 


1. INTRODUCTION 

Spacecraft in geosynchronous orbit exhibit charging when immersed in hot 
substorm plasmas. The potential that a surf-ice attains depends both on the 
material's properties and on the plasma's characteristics. Secondary emission 
of electrons by electrons and ions, as well as backsrattering, determine the 
floating potential. These are both functions of the angle of incidence of particles 
on the surface. 

In this work, the angular distribution of incident particles is obtained for 
attracted and repelled particles for cylindrical spacecraft. Results are given 
for the plane and spherical cases as well. 

For isotropic spacecraft materials, the backscattering coefficients, the 
secondary emission coefficients, and the shielding efficiency depend on impact 
energy and impact zenith angle, but are independent of azimuth angle. According¬ 
ly, only the distribution of impact zenith angles is important. For convenience, 
only protons will be considered, the extension to other charged particles being 
obvious. First the problem will be precisely defined, then a simplified model of 
physical reality set up, and, finally, the distribution of impact angles derived. 
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2. DK.KIMTION Oh PROBLEM 


Kind the flux, j of protons striking an element of area AA with impact 
kinetic energies within the band K t Jlv'2 and with impact zenith tingles within 
the band 0 t A??/2. 


3. MODEL 


The model is that of an infinitivelv long cylinder of radius X . Revotid a 

l ) 

"sheath" of radius R^, the plasma is isotropic and is at zero potential. The 
potential within the sheath is modeled by: 


ln(K/R .) 

o ln(R /R ) 
o s 


( 1 ) 


1'he subscripts "o" and "s" will be used to indicate values of variables at R and 

R .. No subscript will be used for general values of R. 

& -■> - 1 - 1-1 
In the model, the differential flux f(K, protons cm “ see ster eV is 

a function of kinetic and potential energies oniy — except in the "shadow" of the 

Satellite where it is zero. Both potential ($) and kinetic (K) energy are expressed 

in electron volts. Any impact angle and energy corresponding to a trajectory 

originating elsewhere on the satellite is considered to be "shadowed.” 


4. DERIVATION OK IMPACT ZENITH ANGLE DISTRIBUTION 

The impact proton flux is the product of four factors, namely: energy band 
width, unshadowed solid angle, projected area, and differential flux. Tins may 
be w ritten as 

j - (AEKA0 d> sin 0 )<AA cos 0 )[f(E , « )) (2) 

•O 7 . O O O’ O 

where = total unshadowed azimuth angle. The solid angle is that which would 

be subtended at the earth's center by the area enclosed by latitudes D ± At? 2 and 

‘ o 

longitudes ij> ± $ z /2. The azimuth angle will be taken as zero perpendicular to the 
cylinder's axis. If shadowing occurs, it will be at the larger angles. 

Equation (2) may be rewritten as 

j ; (AE A<? AA f(E , 9 )1 [sin 8 cos B ) (4* ) (3) 

J o o o o o ' 


fi 





where b m - maximum unshadowed azimuth angle, <t> being measured So as always 
to lie in the first quadrant. The factor 4 accounts for the four quadrants. A 
further rewriting yields 

i = K <f> sin 2 0 
■o tn o 

K = ^ AO AA f(E , $ ) (4) 

2 o o 

0 < 0 c ff / 2 . 

““ m 

Note: sin 2 0 =2 sin 0 c os 0 
o o o 

Equation (4) is very general. The geometry of the spacecraft's surface and field 
enter into the equation only through $ the maximum unshadowed azimuth angle. 

5. SPECIAL CASES 

For positive spacecraft potentials, the reverse trajectories curve away from 
the spacecraft; therefore, there is no shadowing and 

<1> --- v/2 if 4 ? 0 . (a) 

From Eq. (4), the angular distribution is; 

j = (tr/2) K sin 2 0 . 

■ o ' 

The other special case is that of looking for protons of impossibly low energy. 
Here we may write 

<j> - 0 if 9 < -E , ( ( >) 

m o o 

which has meaning only for negatively charged spacecraft. 


6 . MAXIMUM UNSHADOWED AZIMUTH ANGLE 

The simplest cases have been dealt with in the preceding section. This 
section will consider cases defined by: 





( 7 ) 


-E < $> <0 

o o 



The model has zero axial field. Therefore, the kinetic energy associated with 
the axial velocity component is independent of radial position and may be written: 

E - E sin 2 <j> sin 2 6 , (8) 

a o o o 

where E is the impact energy and <4 o and 0^ are the impact angles. This rela¬ 
tionship follows from the breakdown of the impact velocity into four components, 
as shown in Figure 1. The model has a radial field about the cylinder's axis. 
Thus the tangential (l to radial and axial directions) kinetic energy is subject to 
conservation of angular momentum and is given by: 

E (R) = E (R /R) 2 cos 2 <4 sin 2 6 . (9) 

t o o o o 

The radial kinetic energy is simply: 


E p = E cos" 0 . (10) 

The total proton energy is independent of the radius, and beyond the sheath it 
equals kinetic energy. The total energy is thus: 

E = E+ E+E f + $ , < 0 . (11) 

rat 

At the sheath's edge Eq. (11) becomes: 

E = E cos 2 0 + E sin 2 <4 sin 2 0 + (R /R ) 2 E cos 2 <4 sin 2 0 (12) 

s so o o o' S o o o 

and at the satellite's surface it becomes 

E = E cos 2 0 + E sin 2 $ sin 2 0 + E cos 2 <4 sin 2 0 + $ (13a! 

o oo o oo o oo 


E = E + $ 
s o 


$ < 0 . 


The limiting case corresponds to a vanishing radial energy at the sheath's edge or 


cos 0=0 
s 


I 
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Equations (4) and (lt>) together give the zenith angle (0) distribution of protons 
striking the surface of a cylindrical satellite. 

The angular distribution derivation given has considered only the potential at 
the inner (R - R ) and outer (R = R g ) edges, not the potential within the sheath 
(R < R ,). It may be shown that the potential within the sheath is not controlling, 

s ' i 

provided that at all points within the sheath it is equal to or more negative than 
that of the inverse square potential function: 


$ 


(R /R) 2 -1 
s _ 

(r,/r r- i 




<t < 0 

o 


(17) 


This function is adjusted to match the model at the sheath's edge. The logarith¬ 
mic function of the model satisfies this condition. The proof is similar to that 
given for the spherical case for objects in a gravitational well. The proof is only 
valid if the potential is a function of radial position alone. 

The distribution of impact zenith angles according to Eqs. (4) and (lii) is 
shown in Figure 2 for several cases. The deviations from the symmetrical 
sin (2 0) distribution are due entirely to shadowing. The shadowing is equivalent 
to a "horizon" at less than 90° below the zenith due to refraction of the protons. 



0° 10° 20° 30° 40° 50° 60° 70° 80° 90° 


ZENITH ANGLE 

Figure 2. Distribution of Zenith Impact Angles for Protons on a Cylindrical 
Satellite: (a) Thin-Sheath C ase 
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T. DISCI SS|()N 


I'hi' •hf.iy given here is exact in i licit it is i . • mph ely n msis' • oh the 
model. The model, hinvrver, represents a physical impossibility in that tli« *-!•*» - 
• ric field is discontinuous tit flic surface of the sheath and an approximation ;n 
that end effects are taken to be vanishingly smalt. Only if the nearest end is tit a 
distance very large compared to the radius will its effect become vanishingly 
small. file theory does, however, represent a limiting case and, as such, is 
useful. 


8. SIMM VIO 


The angular distribution of impacting repelled particles is sin (2 '>), as given 
by Kq. (4), since there is no shadowing for repelled particles, electrons inci¬ 
dent on a negatively charged surface and ions on a positively charged surface w ill 
impact in this way. In order to calculate surface potentials, taking account "f 
sec ondary electron emission, one must know the distribution of impact angles a., 
well as the secondary electron yield as a function of tingle. 

The angular distribution of impacting attracted particles (electrons <>n a 
posif’ve surface, and positive ions on a negative surface) is more complicated. 

The impacting angular distribution fin- a cylindrical surface for attracted parti. 1« s 
has been derived here, as shown in Figures 2 (a), 2 (b), and 2 (<•), for a range- of 
energies and sheath thickness. 

The distribution of impact angles (attracted particles) for a plane surface is 

Mi* same as 'he case of zero sheath thickness. For this rase the distribution is 

proportional to sin (2 0) for 0 less than 0 , where 

m 

- cos' 1 (4>/K) 1/2 , 

and where K is the impact energy. If 6 is greater than 0 the distribution func¬ 
tion is zero. For spherical surfaces there is no shadowing, so the distribution is 
again proportional to sin 2 0 in the thick sheath limit. 


12 







